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ROMA  -  1STITUTO  P OL1 071 A FtOO  DELLO  STATO 


SUPERCRITICAL 

GASEOUS  DISCHARGE  WITH  HIGH  FREQUENCY  OSCILLATIONS 

LUIGI  CROCCO  (*) 


K I.4SSUNTO,  -  Ari  conoscenza  del  com/H>rtamento  Ui  un  effusore  aupercritico  quando  lc  eondizioni  a  monte  </♦/• 
V  effusore  sono  osciilatorie  e  importante  in  couni  ssiotie  con  vari  probl  tni  rclativi  at  }>roputson  a  reazione . 
Tsien  ha  rcrerUeiru  nte  analizzato  (a  quest  tone  net  caso  di  piccola  ampiezza  dt  osciUazione  e  (ter  uua  distri- 
buzione  iineare  di  velocitd  nelia  parte  sub  son  tea  delTsffusore.  Egli  ha  dato  la  soluzione  nei  due  cast  cstremi 
di  frequenze  motto  basse  o  >nolto  altc  facendo  uso  di  alcune  ipotesi  rrstrUtive.  Qui  lo  studio  % dene  esteso  ad 
ipotesi  piii  gene rait;  inoltre  net  caso  di  osctllazioni  isentropiche  viene  data  la  soluzione  per  tutte  le  frequenze. 


Introduction 

Many  gus  How  systems,  with  or  without  combustion, 
are  terminated  by  a  discharge  nozzle.  When  the 
stability  of  these  systems  is  studied,  it  is  necessary 
to  determine  the  behavior  of  the  nozzle  under  oscilla¬ 
tory  conditions,  which  may  be  quite  different  from 
the  corresponding  behavior  in  steady  state.  The 
problem  has  recently  been  treated  by  T.sien  (ref.  1) 
for  supercritical  discharge,  this  case  being  the  most 
interesting  for  its  applications  to  the  stability  of 
combustion  in  rockets  or  in  combustors  for  jet  de¬ 
vises.  Tsien  has  investigated  the  case  in  which 
the  oscillations  in  the  incoming  flow  are  isothermal, 
and  therefore  non  ieenfcropic,  and  has  computed  the 
departures  from  the  steady  state  behavior  in  the 
range  of  low  frequencies,  as  well  as  the  asymptotic 
response  to  very  high  frequencies.  Unfortunately 
the  frequencies  appearing  in  the  study  of  high  fre¬ 
quency  instability  are  likely  to  be  in  the  interme¬ 
diate  range,  where  none  of  Tsien’s  solutions  can  be 
used.  Also  the  isothermal  condition  is  not  represent¬ 
ative  of  most  of  the  actual  cases.  The  purpose  of 
the  present  paper  is  the  extension  of  Tsien’s  treat¬ 
ment  to  the  non -isothermal  case,  and  especially  the 
determination  of  the  nozzle  behavior  in  the  interme¬ 
diate  range  of  frequencies. 


Thu  aquation* 


Calling  p,  p,  u  the  pressure,  density  and  velocity 
tn  steady  state,  completely  determined  by  the  shape 
of  the  nozzle ,  and  p  +  p',  p  +  p',  u  +  u'  the  cor¬ 
responding  values  in  unsteady  conditions;  and 
assuming  the  perturbations  p',  p\  u'  to  be  small 
compared  with  the  unperturbed  quantities,  Tsien 
has  written  the  continuity  and  momentum  equa¬ 
tions  in  the  following  form,  retaining  only  the  first- 
order  terms  in  the  perturbations 


\  u  '  '  p  u  I  dx  ^x\  u  i 

^  p’  du  _p_  ?  /P'V 
~  p  dx  pu  1  x\  p  I 


(21 


(*)  Prof— ore  di  Motor!  AerosmUd  alia  Scooted 'r«ecnerte 
Aero  nestle*  dsli'Cnl  verett*  dl  Rom*.  Attuaimente  Robert  H. 
Qoddsrd,  Prof*— or,  Jet  Propulsion  Center,  Princeton,  N.  T* 


x  being  the  distance  along  the  nozzlo  ami  t  the  time. 
The  third  equation  between  the  dependent  va¬ 
riables  //,/>,  p'  p  and  u'  u  is  the  energy  equation  or. 
more  simply,  the  equation  expressing  the  eonataneo  of 
entropy  of  any  fluid  mass  when  wo  follow  its  motion  : 


where  «S'  is  the  entropy  perturbation,  c,  is  the  constant 
volume  specific  heat,  and  y  the  adiabatic  index. 

In  these  equations  u,  duidx  and  p/pu  are  to  be 
considered  known  functions  of  x,  determined  by  the 
nozzle  shape.  Duo  to  the  linearity  of  equations  fl], 
[2],  [3}  the  harmonic  form  of  oscillatory  time  depend¬ 
ence  can  be  chosen,  and,  using  the  complex  repre¬ 
sentation,  the  dependent  variables  can  be  written  as 


*  9  (x)  e1 1 


p 

£  =  8(1)."' 


[<] 


u 


—  v  (x)  e’ 


where  o  is  the  angular  frequency  and  9,  8,  v  are 
complex  functions  of  x  alone.  At  the  nozzle  entrance, 
x  =>  x„  the  three  functions  have  certain  values 
9m  In  w*  The  problem  in  which  we  are  interested  is 
finding  the  distributions  of  9.  8.  v  along  the  nozzle, 
but  expeciaily  determining  the  relations  between 
9,,  8#,  v  These  relations  will  in  fact  constitute  the 
boundary  conditions  to  be  applied  to  the  rest  of  the 
flow  system  as  a  result  of  the  nozzle  presence  (1). 

Equation  [3]  is  immediately  integrated  as 


the  arbitrary  function  /  being  in  general  determined 
by  the  known  time  dependence  of  the  entropy  at 
x*=*  x».  For  the  exponential  time  dependence  assumed 


<H  Instead  of  determining  the—  boundary  condition*, 
Tsien  baa  preferred  to  Introduce  *  •  Transfer  (unction  •  relat¬ 
ing  the  fractional  variation  of  ma—  flow  to  the  fractional  va¬ 
riation  of  pree sore.  The  definition  of  Tsien1*  transfer  function 
can  be  easily  extended  to  the  non  isothermal  case.  However 
we  prefer  not  to  use  this  concept  which  baa  no  direct  utility 
in  the  applications.  The  present  procedure  is  closely  related 
to  the  conventional  use  of  tbo  •  Impadenoe  -  in  acoustic 
treatments. 
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in  equations  [4]  wo  obtain 

V  (•*)  -  Y  S  (*)  =  O  •■*!>  (  —  '  <■'  j ' 


■  i  -  -3)«-(y 

,  I) 


I )  ( I  —  i)  ' 
<@(Y-  II  I 


,/* 


1 1 1 1 


where  the  constant  o  represents  the  amplitude  of 
♦  he  entropy  oscillation  divided  by  »%..  With  the 
assumption  [4]  and  the  relation  [.»)  equations  (!] 
and  [2}  are  rodueed  to  tho  following  system  of  ordi* 
narv  differential  equations  in  v  ami  $ 


/  v  ri 

*3 - r  a  . 

dx  d , 


and  cliiniimt ing  v  the  help  « >t  {111  from  one  of 

the  equations  |t>]  we  obtain 


z{l-z) 


•r-  s  _  / 

dz1  \ 


- 1  ?_  \  <l*  *  * ['2  V  -a)  s 

V  ■+■  I  /  <iz  i  (  v  -y  I) 


id 

(*r* 

\x) 

i  (y  —  1)  +Y  2 

(12) 


'lv  c*  dl  t  dll  i  . .  '/« 

,  , - - - M  2-J - -  <6>  V  —  (V—  I)  —  6  = 

</x  «  d*  I  </x  /  //«c 

=  7^) exp 


[6] 


in  which  c  is  the  sound  velocity.  It  is  easily  seen 
that  these  equations  present  a  singularity  at  n  =  c, 
that  is  at  the  sonic  throat;  where  therefore  only  one 
family  of  solutions  remains  regular. 

Equations  [6)  could  be  discussed  ami  solved  nume- 
rically  for  general  nozzle  shape.  However  the  nu* 
msrical  integration  can  be  avoided  if,  following  Tsien, 
we  confine  our  attention  to  nozzles  in  which  u  increa¬ 
ses  linearly  with  x  in  the  subsonic  portion  of  the 
nozzle,  this  condition  being  not  too  restrictive  since 
many  actual  nozzles  have  practically  linear  velocity 
distribution  near  the  sonic  throat;  and  since  this  is 
the  region  in  which  An  analytical  solution  is  particul¬ 
arly  useful.  Wo  take  therefore 


du  u  cm  rm  ■ —  u, 

dx  x  x+  l 


with  c*  representing  the  so-called  critical  sound  speed 
readied  at  the  throat,  where  .c  ~  x0,  and  l  --  -r#  —  x, 
representing  thj  length  of  tin  subsonic  portion  of  the 
nozzle. 

Moreover  we  take  with  Tsien  a  new  mdtqiendent 
variable 


df-m 


in  terms  of  which  we  have 


cJ 


[9] 


where  u  represents  the  assigned  value  of  [8]  at  the 
nozzle  entrance;  and  we  define  a  reduced  angular 
frequency 


Xm  u  l  ta 
~c*  Cm  —  ut 


00] 


Equation  (I2|  is  a  non  homogeneous  complex 
hypergeometriu  equauou.  with  singularities  at  z  =  0, 
•  «  l,  and  i  —  X-  Of  these  singularities  only  the  one 
at  z  =■=  l  is  important  for  our  problem,  since  the  other 
are  out  of  the  range  of  variability  of  z,  which  musts 
be  contained  between  a  non-vanishing  minimum 
at  the  entrance  of  the  nozzle  and  a  finite  maximum 
at  its  exit.  'The  singular  point  z  —  1  is  a  particular 
case  of  the  more  general  singularity  of  equations  [fij. 


Tha  condition  at  tho  sonic  throat 

Tsien  has  observed  that  since  the  motion  is  super¬ 
sonic  in  the  diverging  part  of  the  nozzle  no  wave 
can  be  transmitted  backwards  to  the  throat,  and 
therefore  the  propagation  of  tho  oscillations  must 
always  take  place  toward  the  downstream  direction. 
He  has  therefore  used  the  condition  that  the  propa¬ 
gation  velocity  (-n  must  always  be  positive.  If  f 
is  defined  hv  the  formula 

;,(;><-£(  f)-” 

substituting  from  (41  wo  have  C  -  -  *  w $l(d&(dx). 
Tho  c  m  lition  that  C  must  remain  of  the  same  sign 
throughout  tho  nozzle  means  that  it  mtist  never 
vanish,  and  therefore  (d  must  remain  finite. 

This  condition  is  not  satisfied  if  the  solution  is  sin¬ 
gular*  with  tho  important  consequence  that  only  a 
solution  which  is  regular  is  compatible  with  the  condi¬ 
tion  of  downstream  propagation.  Hus  conclusion, 
which  is  true  for  tho  most  general  cnee  considered  in 
equations  (fij,  allows  a  more  concise  anil  mathemat  i- 
cally  more  dofinite  expression  of  the  condition  at 
the  throat.  The  same  result  is  obtained  on  a  more 
physical  basis  by  considering  that  a  wave  of  finite 
amplitude  at  the  sonic  throat  cannot  send  but  infi¬ 
nitesimal  waves  upstream,  since  the  upstream  pro¬ 
pagation  velocity  is  zero,  and  therefore  only  waves 
of  infinite  amplitude  at  the  sonic  throat  can  send 
finite  waves  upstream.  The  absence  of  upwards 
moving  waves  is  therefore  connected  with  the  sup¬ 
pression  of  all  singularities  at  the  sonic  throat  (I). 


Introducing  (7J.  [8),  (9J.  (10}  in  equations  [fi]  we 
find,  eliminating  dv  dx: 


(l)  The  Author  U  Indebted  to  dr.  b.-Lesatora  n^r* 

» dtsouMion  of  the  oonditlon  of  noo-ein«ule*ity,  baaed  on  oon* 
LioMt  methods  of  aon-ateady  one-dtmeoatounJ  gaadynamwa. 
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Solution  for  low  fraquaneiat 
In  tha  non-ltotharmal  caio 

Thmmi  Mas  Insiusi  the  ease  of  isnthemml  oscitla* 
tions.  1 1  in  treatment  is  iintneiliately  extended  to 
the  non-isotliermal  ease.  Assuming  a  proscrilteil 
temperature  oscillation  at  the  entrance  of  the  nozzle 


lit'*1 


[13] 


the  expression 


d  »<‘> 
i  it 


-  (log  z  —  :)  — 


e  r  1  . 

T  |  Ylr 


J..  ;_il|+C 


[22] 


whore  the  integration  constant  C  has  to  l»e  dotertnin- 
♦'<1  in  such  n  wav  that  the  right  hand  mile  of  equa¬ 
tion  {22}  vanishes  at  z  I,  ho  that  d  may 

remain  Jlnito  at  thin  point.  We  obtain 


wo  have  between  9,  and  5,  thn  relation 

9,  —  8,  =*  0. 

On  the  other  hand  equation  [ft}  given 

9*  —  y  w  c 

and  therefore  we  obtain  the  relation 


(»*! 


[15) 


[16] 


which  determines  the  amplitude  of  the  density  oscilla¬ 
tion  at  the  entrance  of  the  nozzle  when  0  and  o 
aro  assignod.  Of  course  one  could  prescribe,  instead 
of  arbitrary  values  of  0  and  a,  arbitrary  values  of  9 
and  8,  and  determine  the  corresponding  values  of  6 
and  o. 

In  Tsien’s  case  0*0. 

At  low  frequencies  we  can  expand  all  quantities 
in  powers  of  i $.  since  equations  [II]  and  [I2J  contain 
only  this  combination  of  1  and  fh 

Wo  lake  therefore 

8  <*,*)«*•»<*)  +  «‘3*,,M0+  ....  [17] 

v{*.  3)  *  *♦»(?)  4-  4-  ...  fl«] 


lieplacing  [17]  in  equation  (12)  and  equating  the 
coefficients  of  the  same  powers  of  i  0,  the  equation 
breaks  up  into 


•l*  8<» 

dx* 


#  d  8* 

-*~a r 


a  0 


[!») 


*<l  -«) 


d  J'" 


<1  8* 


[20] 


Tha  aolution  »f  f 1 0|  which  is  non  singular at:  >•  I 
ia  8<*>  c  const.  and  thoreforo.  ainco  tha  aolution  (IT] 
muat  linlil  at  3  ~  II.  from  [ID]  »«  obtain 


8*> 


[21] 


Introducin'*  thla  value  In  equation  [20],  ami  Inte¬ 
grating  thla  first  nrilat  aquation  in  d  *i"/ifx  wo  obtain 


d  81b  0  log  :  4*  I  —  * 

=  YJ  -  I  ( I  —  .„3. 

_ o  flog :  +  1  —  : _  T  — t|  1 

-  Y  (Y  —  1  >  I  <  1  —  «)*  *  ]’ 

Tt  is  immediately  checked  that  this  expression  in 
regular  nt:  a  1.  An  additional  integration,  with 
the  condition  following  from  (16)  and  [21] 

8<i>,  *0  (24) 

would  give  the  (non-singular)  expreaaion  for  8<*>  (*). 

Replacing  now  [17],  [18]  in  equation  [11]  and  Again 
equating  tho  coefficients  of  the  same  powers  of  ifl 
one  obtains,  with  the  use  of  (21),  (23),  (24) 

v«*>  n  V*1  13  -y  ; 


V— 4^- ■!?!*■)♦ 


and  therefore,  recalling  (17),  (18): 


v#  _  v/**  4-  i  (J  v/1*  4*  ...  */  —  I 

8.  =  8/*r+TpT8#.,,'-u  ...  *  T 

+  t=±. 

9  —  o 


i »  +  ....  (3  <T) 


[2#] 


[29] 


This  quantity,  representing  tho  complex  ratio 
Iwtweeu  tho  frnctiunal  variationa  of  velocity  anil 
density  at  the  entrance  of  the  nozzle,  and  analogous 
to  the  reciprocal  impedance  of  acoustic*,  can  Be  used 
as  Boundary  condition  for  the  root  of  the  flow  system 
in  oscillatory  state.  Wo  see  that  [2fl]  depend*  only 
on  3.  «  anil  tho  ratio  fl/a.  The  flrst  term  of  tha  aeriee 
(29)  applies  to  the  steady  state  and  contains  as 
particular  cases  the  ieothermsl  cose  9  —  0  where 
the  velocity  cannot  rhango  (w  —  0):  the  ixentropio 
case  e  —  0,  w/8,  —  (y  —  l)/2.  and  the  ieopicnie  case 
9  ~  a.  8,  «•  0.  Equation  (28)  ehows  that  even  for 
modorate  9  the  boundary  condition  impoeed  by  the 
nettle  can  ohange  considerably  sine*  both  the  phase 
and  the  amplitude  ratio  between  velocity  and  density 
fluctuations  are  alTected  considerably. 

The  same  procedure  can  be  used,  without  sub¬ 
stantial  difficulty,  to  compute  higher  order  terms 
ir  i  0. 
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Solution  for  high  fraquoncioi 
In  tho  non-l*otharmal  case 

Following  n  procedure  similar  1"  llio  one  used  by 
Tsict  «•«  first  dotertnino  tho  particular  solution  S* 
of  equation  [1-1  taking 


*•“>-* w(t)  ’ 


where  the  constant  C  is  determined  by  tho  condi¬ 
tion  [HI]  at  z  =•  z,  and  is  oxpnwHutl  by 


C  0-7)+- 


and  replacing  it  in  equation  f  1  -!•  An  equation  for 
7.  (:)  is  obtained  whicli  can  he  solved  by  series  taking 


Z  (r)  =.  Z'»  (z)  +  y-  2<> I  (:)  +  .... 


and  equating  tho  coefficients  of  the  same  powers  of  i  p. 
The  result  is 

— ■ -7(iP(,  +  tfiTrp5!  +  -)-OT 

Replacing  [27]  in  equation  [11],  the  correspond¬ 
ing  value  of  V  is. round  to  be  zero  (by  pushing  the 
expansions  to  higher  powers  it  is  actually  found 
that  tho  first  non  vanishing  term  in  tho  series  for  v 
corresponding  to  the  particular  solution  [27]  is  the 
term  in  (*  ?)-*).  The  *•  fulfils  the  condition  of 
regularity  at  z  =  1,  hut  not  the  condition  [Iff] 
at  z  —  r,. 

Therefore  a  solution  of  tho  homogeneous  equa¬ 
tion  corresponding  to  [12]  must  lie  determined. 
If  wo  take  8  exp  (»'  P  X  <:))  and  replaco  in  this 
homogeneous  oquatinn,  wo  find  that  tho  derivative 
,1  XI, It  =*  y  (Z|  satisfies  the  lticcati's  equation 

Introducing  in  this  equation  the  series 

y  (z)  «■  y'*1  (z)  +  ^y"'(zM-.... 

and  equating  the  term*  with  equal  power*  of  we 
obtain  ^ 

„..i  „ _ L - (|  'l  +  X±1  _  l)  [28] 

V  (y-H  l){»  — *> VI  %  z  / 

r,  Y-t-l  -|„_ag- 

r  r+i  l  —  l  ,lr 

Sr1  +  _T““  t«j 

_Sv,« _ 1 _ LI 

r  Y+  I  «  J 

where  only  the  solutions  remaining  regular  at  t  -  l 
have  been  considered.  The  solution  of  equation  [12] 
can  be  put  now  under  the  form 

ft  »«0exp|  ifij9  -f  +  4-  ft*  [30] 


Coming  now  to  oquntion  [11]  and  recalling  tlmt  S* 
does  not  bring  any  contribution  to  v  up  to  terms 
of  order  (i fj)—*,  up  to  this  order  wo  can  express  v  a* 

v  =  + 7Lr.,,,  +  ".) 

Coxp(rpJ  y'^tlz  -f  /  ytv  dz  +  -j  [32] 

and.  after  substitution  in  equation  [II]  and  compa¬ 
rison  of  terms  with  same  powers  of  » {1,  wo  find 

,».|,1+y+iL|jr  r3j] 

1)“'  -  Y—  1  —  2n'"  —  <Y+  1)(1— *)»'n 

y'“  boiug  givon  by  [28],  [20]. 

Finally  computing  v  at  z  »»  z,  from  [32]  and  recal¬ 
ling  tho  condition  [10]  wo  find 


»•  „  _ Y  1  /Y-l 

8.  •  -a*  f  ip  I  Y  »“ 

0--  \ 

i  +  -  (3 >  I) 


i  /y— 1  _* _ n-'*1  . 

p  (  y  8  — oy  +  J  *  Y~l} 

\  ~T.  2 


with  Tj.'11  given  by  equations  [33],  [28],  [29]  at 

Higher  order  terms  in  expansions  [27],  [30],  [32] 
and  [34]  can  easily  bo  cimn'Oted.  Again  wo  see 
that  v,/8,  do|ienda  only  on  p.  -  and  6  o-  I'or  v  0 
and  [J  *  70  [34]  gives  Tsi,  n'a  -exult. 


Practical  i 


It  is  interesting  to  investigate  if  equation*  [26J 
or  [34]  can  bo  used  for  stability  computations  of 
practical  combustion  systems.  It  seems  to  day  esta¬ 
blished  that  high  frequencies  ap|iearing  in  unstable 
combustion  are  close  to  the  natural  modes  of  oscilla¬ 
tion  of  the  combustion  system,  that  is  —  for  in¬ 
stance  —  to  the  organ-pipe  frequencies  if  the  com¬ 
bustion  chamber  is  of  elongated  form.  Tho  funda¬ 
mental  organ-pipe  mode  has  an  angular  frequency 
o  —  n  c ,/£,  e,  representing  the  sound  velocity  in 
the  ohamberand  L  its  lenght.  Hence  from  [10]  we 
have  p  -  ir  e,  J/(%  —  u,|  l.  Xow  c,  (c.  —  «.)  is 
clous  to  one,  so  that,  if  the  length  of  the  subsonic 
portion  of  the  nozzle  is  around  1/1  of  ths  length  of 
the  chamber,  p  Is  around  unity  for  ths  fundamental 
mode,  around  two  for  the  second  mode  end  so  on- 


L’ABUOTECMICA 


|V..I,  XXXIII.  las.-.  l\  1  | 


This  approximate  computation  shows  therefore  that 
unless  tho  geometry  of  the  system  is  quite  different 
from  the  assumed  one  tho  resists  of  the  previous 
seel  ...ns  lilt  not.  Kite  imieli  help  if  the  fre.|Ueiu-ies 
are  eloso  t..  the  naturui  imsles  of  psrillntion  of  llio 
elmmls-r,  as  —  lor  instance  —  in  the  study  of  *  Korean. ■ 
inir  ■  eou.litions  in  ro<*ket  motors. 


Solution*  for  all  froquanela* 

In  tho  itontroplc  caio 

Fortunately  it  is  (sissiblo  to  tako  advantage  of  the 
analytical  pro)K>rties  of  equation  [12]  to  extend  the 
eoniputatiojUL  tn-tlis  whole  range  of  frequencies!  Tu 
order  Vo  avoid  tho  oomplicatitmconnoctod  with  the 
determination  i>T*tlio  particular  sonimm  ol  equa¬ 
tion  1 1 21.  we  have  hero  considered  only  me  laontrojitu 
csae!~^rMch  proEitity'~repre»eiiru  quite  cloSiTytho 
actual  t.Tienoinenon.  W  ith  a  »  o  equation  Q~-l  IH 
reduced  to  an  homogeneous  hvpertreoinetrie  equa- 
tiun.  of  which  the  solution  remainiuir  non-singular 
at  i  -  I  is  Riven  l>y  the  known  hypergcoinetrie  series 
in  powers  of  1  — •  z: 

F{n,b\c;  X  —  ^  1  -r  —  (l  — *)*f 

.  r»<«»—  1)6(6-+-  1>(I  —  *)* 

+  e  (c  +  I) - 5T“T-  l30] 

with  n.  b.  c  given  by 

e.»  +  HI»l[l  + 

'  Y+  [Ml 

„6_-ii/,  +  L_Y 


In  principle,  then-fore,  mir  prohli-in  is  solved  taking 

t  —  C  F  (o.  6;  c;  I  —  r)  [371 

and  eoinpiitimt  v  from  equation  ( 1 1 1  with  e-"0  and 
with  d  F:dz  nhtaine.1  from  the  differentiation  or  [35], 
Tho  ratio  v,/8.  is  therefore  dcterminoil  and  is  itlile- 
pendent  of  C. 

However  thie  procedure  rannot  bo  followed  in 
practical  ease*  sineo  z,  is  generally  quite  small  (I) 
and  I  —  z,  is  close  to  one,  in  a  region  where  the 
convergence  of  [3.1]  is  too  poor.  Tlie  difficulty  can 
lie  overcome  with  tlie  help  of  the  properties  of  the 
solutions  of  the  liyporgeometrie  equation*.  .Simpler 
<levelo|Mnenta  are  obtained  using  Instead  or  series  [38], 
the  pniportlonal  series 

/  (rs,6iCi  I  —  t)  r  tpFJ ■■  f  «».  *><  et  I  —  s)  — 

.  .. 

“,?,<l“,^P|.+  i)Ti:+-.i 

where  T  (r|  represents  ihe  gamma  function  of  argu- 

tu  It  ttw  Maek  nasibtf  st  Ike  eettl*  eettaees  Is  .(,  I  s, 
Is  erased  .M.  Eras  (nr  Msefc  ssstSs*  I-*.  Is  skw*  to  .ft. 


moot  x.  This  particular  solution  can  be  expressed 
as  a  linear  combination  of  tlie  series  expressing  two 
particular  solutions  of  tho  liypct-gcometric  equation 
in  |M>wcrs  of  r. 

One  such  relationship,  particularly  useful  in  tho 
present  case,  is  .ref. 

/  (o,  6;  ci  l  —  z)  = 

„_._i  ;/|r  —  r  —  h:  I  +  £  —  >»  —  *)_ 

I'(r  —  .<)  I'  (c  —  b)  cos  (c  — «  —  b)t r 

[30] 

where 

t/{a\b'ic';z)  =s  —  it  cot  |c'  “1  [/  (o',  6':  c':  z)  — 

-:‘—'fln’+  I— c',6'+l-e';2— e':t>] 

1 4H] 

the  two  functions  /  Is-inc  given  by  the  cortes|Mind- 
ing  series  |3H].  With  the  particular  valuo  13(1]  of  r, 
the  value  of  c'  to  use  in  [40]  is  given,  following  [311], 
by  3.  Vow  for  integer  c'  «  n  the  quantity 
in  tho  brackets  of  equation  [40]  vanishes  (as  it  is 
ensily  checked)  and  the  factor  preceding  the  braekotB 
becomes  infinity.  Tho  corresponding  valuo  of  g 
can  be  found  through  a  limiting  procosa  to  bo  (ref.  2), 
for  «  >  I , 

0  (o',  6';  n: z)  a  —  I  («’,  l/i  »:  z)  log  s  + 

rt—  .iris+u'iris-bto 

,r(s^,oris  +  tj)[i))(  ,)  + 

,r.  r(»+i)r(*+«i iyi 

-f-$l*-*-6')— $(»+  I)— +(« +•")]  [41] 

where  <[i  if')  re|iresents.  as  usual,  the  logantlimio 
derivative  of  ihe  gamma  function  with  n-sfss-t  to 
the  argument  r. 

Finally  from  equations  |3X],  [3!l|,  1 4 1]  with  o,  b.  n 
given  by  [.111]  and  therefore  u  —  -  wo  find 


f  (a,  b:  r:  1  —  z)  «  log  z  ^  .-I,;1 + 


-  i]  .1,  />.=■  + 


2(y  +  l> 
i  [1  (2  +  i  >) 


with  tho  cootlioioiits  .1,  given  by  tlio  recurrence 
formula 


I- 


2(y +!)#(«-(- I)  <Y+  R" 


f),-i(i(*  +  .i)  +  +  (*  +  »)  — +  (»)—']'(«+  ')•  [43] 

Tito  computation  of  i),  can  be  performed  with 
the  help  of  the  series  expression  of  i]>  (*): 

+  (*|  —  log  *  —  —  +  —  >  — 7#f“‘ 

r  m  I 

where  B,  -  I  /(I.  B,  -  l/JO .  are  the  Bernoulli's 

numbers. 
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We  obtain 

A,  2  [  s  a  -r  I  *  b 
*  H* 

—  —  (*+•  l)“  *']■ 


±i± 

L 

d  ^+.J 

*f  i'.y  —  6)_,r  — 


The  convergence  of  the  lost  series  is  very  fust  for 
high  values  of  a.  O,  has  been  therefore  computed 
from  [44]  only  for  the  highest  value  of  a  needed  in 
the  evaluation  of  [42j;  and  for  the  other  values  of 
f),  from  the  recurrence  formula 


V.+  ,  -X>.  =  — •  H - l-y - ~  - 

n  4-  a  a  -j-  0  * 


The  solution  [37]  can  be  expressed  as 


8  b:c;l  —  z) 


differentiation.  The  real  and  the  imaginary  part 
of  this  quantity,  the  first  divided  by  and  the 
second  by  [i  for  convenience  of  scale,  ore  given  for 
Y*®  1.2  in  figures  1  and  2.  The  number  of  terms 
used  in  the  computation  are  sufficient  ro  give  very 
accurate  values  tip  to  s  =  .2  for  senes  [42]  and 
reasonable  accurate  data  up  to  :  .3:  the  corre¬ 

sponding  limits  for  series  [33]  are  z  =*  . S  and  2  =  .7. 
I  lot  ween  .3  and  .7  the  dotted  curve  is  only  interpo¬ 
lated.  The  lines  for  fl  _«  0  are  entirely  computed 
from  the  equations 

(  V  4-  1  V  Z  log  2  (log  2—1  — Z)  j 

-  T  -J-)  »««* - -7”i  —  j 

4-0 1  'i  '  s  dz  )J  (Y-f  1)(1-*)* 


with  the  constants  C'  and  C  connected  through  [38]; 
and  [42]  provides  the  series  expansion  suited  for 
computations  at  z  =  z..  Actually  to  solve  the  pro¬ 
blem  of  the  boundary  condition  equivalent  to  the 
presence  of  the  nozzle  we  need  only  the  calculation 
of  v/8  and  therefore,  from  equation  [11],  of  the 
quantity  <  1//)  ( dfdz )  =  (1  jF)  (rfFfdz),  independent 
of  the  values  of  the  integration  constants  C  or  C\ 
and  easy  to  calculate  from  the  inverse  ratio  of  se¬ 
ries  [42]  or  [35]  and  the  series  obtained  through  their 


which  can  lx>  derived  with  the  procedure  used  pre- 
viously  for  small  p. 

For  the  first  of  the  two  quantities  goes 

to  zero,  but  the  second  one  takes  the  expression 


»/‘0)  being  given  by  equation  [28].  Both  quantities 
tend  logarithmically  to  infinity  at  z  =  0,  a  value 
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which  can  never  occur  in  practice.  .The  Iwhavior 
of  (I  8.)  if/ 8  dz)4  is  better  seen  m  liguro  ‘.I  where 
this  quantity  divided  by  —  /  3  is  represented  in  the 
o«in|>lt’X  | •lane  for  various  values  at  c,  =-  lt/r  *\ 

From  < I  5.)  ('/$  ite),  anti  equation  ill)  the  quantity 
v,/8,  has  been  computed.  Figures  4  ami  .“»  show 
the  real  end  the  imaginary  part  of  «:s  a  funetion 
of  p  for  different  vahtes  of  r„.  Figu  e  •>  given  a 
represent  at  mu  of  v,  /£,  in  the  roinjilcx  plane,  lrom 
which  the  pluvs'  anti  amplitude  relation;  hetwt'cn 
the  velocity  a*nl  tiensity  fractional  fluctuations  arc 
immediately  deduced.  For  a  given  Mach  imiidier 
at  the  nozzle  entrenee  the  ratio  of  the  amplitude? 
increases  stoat lilv  with  increasing  frequencies.  while 


3. 

i 

1 

\ 

l 

1 

1 

j 

9 

j 

the  phase  m*«'N  tV.  m  zero  »«»  «  maximum  which  is 
always  an  important  fraction  of  -  2  anti  then  tie* 
creases  back  to  zero. 


'Hie  computet l  values  of  v.,’8,  are  sulHcient  to 
solve  problems  in  the  complete  range  of  practical 
frequencies  anil  Mach  numbers  m  the  isentropic 
problem  anti  for  linear  velocity  i list ri but  ion.  Addi- 
tional  computation  art'  requiretl  in  other  cases. 

The  isentropic  ru**'  with  other  velocity  distribu¬ 
tions  ran  bo  coinput  »d  st  irtmg  from  the  solution  for 
the  linear  distribution,  which  is  still  available  near 
Css  I  provided  dirdxz^i  0.  and  continuing  with  the 
numerical  integration  of  equations  |t?l  with  a  =  o.  In 
particular  at  the  sonic  point  r=  1  wo  have  the  relation 


*•—14-/3 
2  4-  ip 


that  is  the  same  equation  given  *>’k  fill  for  finite 
(1/8)  (d  8/dt)  at  z  «  I.  with  p  =  u/(rfu/rfx ),•»,*.  fn 
terms  of  P  the  relation  between  v  ami  8  at  the  sonic 
point  is  independent  of  the  nozzle  shape. 
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It  in  iutere^timr  to  olmerve  that  tlio  fractional..  and  llw  amplitude  of  the  oscillation*  can  I»o  related 


oscillation  of  Mach  number,  jriven  in  tHo  laontmpic  to~th©  amplitude  of  oacillatian  of  pmmiin*,  density 

rase  by  v  •  (y^~  iTT^TTh  at  the  throat  cijual  to  7»F‘loinpernturo  at  the  rm/.zlo  entruneo. 

Finally  wo  observe-" with  Tsirn  that,  rs|treinlly 
\  ' — M  %  —  y)  •  3  for  short  approaches  to  the  throat  the  onedimen* 

2  j  “  2(2  r  i‘;i)  sionnl  assumption  ran  lie  very  much  in  error,  and 

that  additional  error*  are  in  lie  expected  from  the 
re  different  from  zero  for  Unite  The  fart  that  the  fluid  hits  been  roiisidercd  a*  a  frietiam* 

herefore  oseillntitur  around  the  throat,  lo-w,  in-  atim*.  non  reart  ini:  uas. 


[/v,  Y~ 1|  .  <^-Y)-? 

*•*•’•!*  2  i  — 2(2 

ami  is  therefore  different  from  zero  tor  finite  The 
sonic  line  ta  therefore  itseilluf  ititf  around  the  throat, 
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SUPERCRITICAL  GASEOUS  DISCHARGE  WITH  HIGH  FREQUENCY  OSCILLATIONS 

by  l.  CHOCCO 


l 


For  several  problems  connected  with  jet  propulsion  nozzle.  I’tulrr  certain  n'mrictiw  assumptions  he 

it  is  riH|Uir*i|  to  know  the  behavior  of  a  siipercrit teal  has  given  the  solution  in  the  two  extreme  cases  of 

discharge  nozzle  when  the  condition*  iqisttvain  of  very  low  and  of  very  high  frequency.  Here  the 

the  nozzle  are  oscillating.  Tsien  has  recently  <tudy  is  extended  to  more  general  a  sumptions, 

analyze!  the  question  in  the  case  in  which  the  oaciiia-  Moreover,  in  the  ease  m  which  the  iqwtream  oscil- 

tion  amplitudes  are  small  and  the  velocity  axial  1st  ions  are  jsemropie.  the  solution  in  given  in  the 

distribution  in  linear  in  the  subsonic  portion  of  the  complete  range  of  frequencies. 


